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Abstract 

Using a cell dynamic system (CDS) simulation scheme, we investigate the 
phase-ordering dynamics of non-conserved 0{n) models without topological 
defects, i.e. for n > d + 1 where d is the spatial dimensionality. In partic- 
ular, we consider zero-temperature quenches for d = 2, n = 4, 5, and for 
d = 1, n = 3,4, 5. We find, in agreement with previous simulations us- 
ing fixed-length spins, that dynamical scaling is obtained, with characteristic 
length L{t) = t 1 / 2 . We show that the asymptotic behaviour of the structure- 
factor scaling function g{q) is well fitted by the stretched exponential form 
g(q) ~ exp(— bq s ), with an exponent S that appears to depend on both n and 
d. An analytical treatment of an approximate large-n equation for the pair 
correlation function yields g(q) ~ q _ ( d_1 )/ 2 exp(— bq), with b ~ (Inn) 1 / 2 for 
large n, in agreement with recent simulations of the same equation. 
64.60.Cn, 64.60.My 
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I. INTRODUCTION 



The phase-ordering dynamics of systems quenched from a high-temperature disordered 
state into an ordered state is a problem of great relevance in the description of out-of- 
equilibrium pattern formation ffl]. One well established property is the onset of dynamic 
scaling, where the late-time behaviour of the order-parameter correlation functions is de- 
scribed by scaling forms with a single time-dependent length scale L(t). Thus the real-space 
correlation function is found to have the scaling form 

C(r,t) = f(r/L(t)) , (1.1) 

while its Fourier transform, the structure factor, has the corresponding scaling form 

S(k,t) = [L(t)] d g(kL(t)) . (1.2) 

Conventional experimental systems such as binary alloys and binary liquids are described 
by a scalar order parameter. Recently, however, there has been much interest in systems 
with more complicated order parameters such as n-component vectors (the 0(n) model) 
15 1 and traceless symmetric tensors (nematic liquid crystals) ||16|| . In this paper we restrict 



discussion to the 0(n) model. 

Much numerical and theoretical effort has been devoted to understanding the basic prop- 
erties of systems that can support singular topological defects, i.e. systems with n < d. The 
presence of such defects leads to a generalization of the usual Porod law for the large-g tail 
of the structure-factor scaling function g(q), namely g(q) ~ q~\ d+n ) for g > 1 |. This 
result is geometrical in origin, and is independent of whether or not the order parameter is 
conserved by the dynamics 0. 

Very recently, the cases n = d+1, for which nonsingular topological textures occur, have 
been studied numerically (for d = 2) and analytically (for d = 1). Rutenberg and Bray || 
found that the d = 1 XY model (n = 2) exhibits a scaling violation due to the existence of 
two relevant length scales: the phase coherence length and phase winding length. On the 



2 



other hand, the two-dimensional Heisenberg model (n = 3) with non-conserved dynamics 
also violates dynamic scaling due, it appears, to the existence of as many as three separate 
length scales, related to individual texture size, the typical separation between textures and 
the typical distance between textures of opposite charge ||. These texture systems are, 
perhaps, the most complex of the phase ordering systems. 

By contrast, systems without topological defects (n > d + 1) seem relatively straight- 
forward. There is good evidence for the simple scaling behavior described by ( P-.1D and 
( |1.2|) , with characteristic scale L(t) ~ t 1 ^ 2 for nonconserved dynamics. The energy scaling 
approach of Bray and Rutenberg |17| shows that, provided scaling holds, L{t) ~ t 1 / 2 is 



indeed correct for n > d + 1 systems with nonconserved dynamics, and gives L(t) ~ t 1//4 for 



conserved dynamics, again nicely consistent with simulation results [|T0|] and the Renormal- 
ization Group result of Bray 



Much less is known, however, about the form of the structure-factor tail for n > d + 1. 
The recent simulations of Rao and Chakrabarty (RC) [[Hj, with conserved dynamics, for 
the cases d = 1, n — 3 and d = 2, n = 4 show 'squeezed exponential' behavior [i.e. 
g(q) ~ exp(— bq 5 ) with 5 > 1]. In this paper we concentrate on systems with n > d + 1 
and nonconserved dynamics. We consider the cases d = 2, n = 4, 5 and d= 1, n = 3,4, 5. 
In each case we confirm the expected t 1 ^ 2 growth, and find 'stretched exponential' behavior 
[i.e. g(q) ~ exp(— bq s ) with 5 < 1] for the tail of the structure factor, with an exponent 5 
that appears to depend on n and d. 

In an attempt to understand the origin of this tail behaviour, we present an analytical 
approach based on an approximate equation due to Bray and Humayun (BH) fl2"| , which is 
itself based on the 'gaussian auxiliary field' (GAF) method []14| that describes rather well 



the form of the structure factor for nonconserved systems with singular defects (n < d). 
For those systems, this method reproduces, in particular, the generalized Porod tail. For 
n > d + 1, the physical basis of the method is less clear. However, the simple truncation 
of the equation at leading order in 1/n, proposed by BH in another context ]T2"|| , leads to 



an exponential decay of g{q), modified by a power-law prefactor for d > 1. It is noteworthy 



that the asymptotic behaviour is nonanalytic in 1/n: for n strictly infinite, the gaussian 
form g(q) ~ exp(— 2q 2 ) is obtained. The exponential asymptotics of the BH equation were 
noted in recent numerical studies by Castellano and Zannetti |13[ . 

Using a 'hard-spin' model Newman et al. H studied numerically the dynamics of one- 
dimensional systems without defects for n = 3,4 and 5. Measuring only the real-space 
correlations, they found that dynamic scaling is obeyed with characteristic length L{t) = t 1 ^ 2 . 
Moreover, they found the real-space correlation function was very well fitted by a gaussian 
form, which is the exact result in the limit n — > oo. The Fourier space analysis presented 
here, revealing stretched exponential tails, shows that the good gaussian fits achieved in real 
space are misleading. 

Our main results can be summarized as follows: (i) For all our models the characteristic 
length scale required to collapse the data for the real-space correlation function and structure 
factor, is Lit) = t 1 ^ 2 , in agreement with theoretical predictions fl7fl ; (ii) The asymptotic 
behaviour of the structure factor is well described by a stretched exponential of the form 
g(q) ~ exp(— bq s ), where the exponent 5 apparently depends on both n and d and seems to be 
different from the value obtained for the corresponding system with conserved dynamics ||TU 



(iii) An analytical treatment of the approximate BH equation, expected to be valid at large 
(but finite) n, gives a simple exponential modified by a power, g(q) ~ q-( d ~ 1 )/ 2 exp(— bq), 
with the same asymptotic form for conserved dynamics. 

The rest of the article is organized as follows. In the next section, we introduce the 
CDS model based on the time-dependent Ginzburg-Landau (TDGL) equation for a zero 
temperature quench, and we describe the corresponding numerical procedure employed in 
the simulation. Section 3 presents simulation results for a vector order parameter with 
n = 4 and n = 5 components in one and two dimensions and the one-dimensional 0(3) 
model. For the d = 2 systems, we also present data for the real-space correlation function to 
demonstrate the dynamic scaling behaviour. We then discuss the procedure used to obtain 
the asymptotic functional form of the structure factor tail. Next we compare the data to 
the results of the approximate analytic theory. Finally we make some concluding comments 



on our results and a give a brief summary. 



II. MODEL AND SIMULATIONS 

The dynamic evolution of a non-conserved vector order parameter (model A) with n 
components = (0i, 02, • • • , 4> n ), for a zero-temperature quench, is described by a purely 
dissipative process defined in term of the following TDGL equation: 

dffrt) = 5F($(x,t)) 
dt 5$(x,t) ' 

with T a kinetic coefficient that we will set equal to unity, and F the free energy functional 

which generates the thermodynamic force, 
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F[f(x,t)] = J d d x 
with the potential defined as 



^(V0(x,t)) 2 + ^(0(x,t)) 



(2.2) 



V(0(x,t)) = i(l-0 2 (x,t)) 2 (2.3) 

where 2 = Yh=x 0?( x > 0- The ground states, or fixed points of the dynamics, are deter- 
mined by the condition 2 = 1, which defines a degenerate manifold of states connected by 
rotations. In the internal space of the order parameter, this manifold is the surface of an 
n-dimensional sphere. At late times the order parameter is saturated in length (i.e. lies on 
the ground state manifold everywhere). Then the dynamics is driven by the decrease of the 
free-energy associated with the term (V0) 2 in fl2.2|) , through a reduction in the magnitude 
of the spatial gradients. 

We can construct an explicit numerical scheme for the simulation based on a computa- 
tionally efficient algorithm, namely the cell-dynamic system (CDS) ||19|| , which updates the 
order parameter according to the rule 



0(x,t + l) =H(((>(x.,t)) + TD 

with 



-]T0(x',t)-0(x,t) 



(2.4) 



Htffct)) = 0(x,t) + r0(x,t)(l - (0(x,t) 2 ) , 



(2.5) 



where z is the number of nearest neighbors, and r and D are parameters that we choose to 
be t = 0.2 and D = 0.5 in our simulations. 



in the values of the parameters r and D. The CDS is an Euler-like algorithm and for 
convenience in our analysis of the results we use a unit of time equal to the update time step 
r. It should be noted (see Figures 1 and 3) that the scaling regime is reached very quickly 
in these systems without defects, and very long runs are not necessary. 

The two-dimensional systems consist of a square lattice of size 256 x 256 with periodic 
boundary conditions. The physical quantities are calculated as averages over 20 independent 
distributions of initial conditions. The one-dimensional systems have L = 16384 sites (with 
periodic boundary conditions) and we average 100 independent runs. The initial conditions 
for the order parameter components (pi were randomly chosen from a uniform distribution 
with support on the interval (-0.1,0.1). 

A quantity of interest that is computed during the course of the numerical simulation in 
the two-dimensional models is the two-point real-space correlation function 



where < • • • > stands for the average over the set of independent initial conditions (or 
'runs'). A spherical average over all possible distances r = |r| is performed to find the 
isotropic real-space correlation C(r,t). The other function of interest, calculated for all the 
models, is the structure factor, 



We also make a spherical average over all possible values of k with given k — |k|. 

In the calculation of these quantities at each time, the data are 'hardened' by replacing 
the order parameter at each point by a unit vector in the same direction (the fixed point of 
the CDS iteration being a vector of unit length). This procedure accelerates the entry into 



The above numerical procedure is identical to that used by Toyoki differing only 




(2.6) 




(2.7) 
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the dynamic scaling regime, and helps us to elucidate the proper nature of the asymptotic 
tail in the structure factor. 

III. RESULTS 

Dynamic scaling is observed for all the models studied. The scaling regime is reached 
at quite early times, in agreement with previous studies. We show that dynamic scaling 
holds in the two-dimensional systems (n = 4,5), using the characteristic length L(t) = t 1 / 2 
deduced from theoretical considerations ||17|| . This agrees with earlier simulations of Bray 
and Humayun using 'hard-spin' dynamics |$|. 

Figure 1(a) presents a plot for d — 2, n — 4 of the correlation function ( |1 . 1|) as a function 
of distance r for several times, while in Figure 1(b) we show the collapsed dynamic scaling 
function when the analysis is made using the scaling variable x = r/L(t). As can be seen 
from the Figure, the scaling function f(x) is a monotonically decreasing function with the 
generic featureless shape that is characteristic of non-conserved 0(n) models. 

It is of some interest to investigate the small-x behavior of the real-space scaling function 
f(x). In systems with n < d, the existence of singular topological defects leads to a non- 
analytic term of the form \x\ n (with an additional In a; factor for even n), which leads 
to the fc _ ( d + n ) Porod tail in Fourier space |||7|. in the present case, where n > d, we 
expect no such short- distance singularities. Therefore, we consider an expansion of the form 
f(x) = 1 — ax 2 + /3x A ■ ■ ■. In table 1 we present the parameters a and (3 determined from 
the simulations in the range x < 0.5. The ratio r = f3/a 2 should be a universal number for 
given n and d. It will be seen from table 1 that this ratio has the value r ~ 0.59 for n = 4, 
different from the value 1/2 obtained for a gaussian function, which is the exact result for the 
limit n — > oo. For n = 5, table 1 gives r ~ 0.49, already consistent with the large-n result. 
However, in the absence of any short- distance singularity, the small-x behavior provides no 
useful information on the nature of the tail in the structure factor. Consequently, it is more 
convenient to investigate directly the simulation results of the structure factor and extract 
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from them the asymptotic behaviour. We shall see that the behavior in Fourier space is 
clearly non-gaussian, even for n — 5. 

As expected, given the absence of topological defects, the results indicate (Fig. 2) that 
the decay of the structure factor is clearly faster than a power law, in contrast to the 
interpretation of his own similar results by Toyoki |T3|. In order to demonstrate that the 



tail is well described by the stretched exponential form 

g(q) ~ Aexp(-bq 5 ) , (3.1) 

where q = kL(t) is the scaling variable in momentum space, we attempt to find the corre- 
sponding power S in the exponential by plotting lng versus q s and adjusting the value of 5 
until the best linear behavior is obtained in the regime q > 1. During the fitting procedure 
the other two parameters of the fit, A and b, are readily determined. The criteria used for the 
optimum fitting is based on the Pearson correlation coefficient (PCC), which measure the 
strength of the linear relation among two variables and varies between -1 (perfect negative 
linear relation) and +1 (perfect positive linear relationship). We proceed as follows: first, we 
choose an exponent 5 and then perform linear regression; next we change 5 until the PCC 
reaches its maximum value. The regression coefficients are calculated using the values of 
the scaling structure function at the last two times in the simulation. The optimum values 
for system with n = 4 components are 5 = 0.435 with a Pearson coefficient of —.999998. 
The other two parameters are In A = 13.21 and b = 8.19 This result is presented in the Fig. 
2(b). 

We turn now to the description of the case n = 5, following a similar analysis to the n = 4 
model. Fig. 3(a) shows the correlation function as a function of distance for different times. 
In this model the collapse is also achieved using the characteristic length L(t) = t 1 ^ 2 , as can 
be observed in Fig. 3(b). Therefore, both models are consistent with dynamical exponent 
z = 2. The corresponding scaling plot for the structure factor is shown in Fig. 4(a). A 
more important effect is observed in the structure factor tail, because in this case it has 
also a stretched exponential but with an apparently larger exponent. Following an analysis 
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similar to that used for n = 4, we find that the value of the best fit value of the exponent 
is 5 = 0.613, and the corresponding PCC in the regression is —0.999998. The other two 
parameters are hiv4 = 7.57 and b = 4.39. In Fig. 4(b) we plot In g against q 6 and the linear 
behaviour is clearly seen. 

Comparison between the real-space correlation functions of the n = 4 and n = 5 models 
shows that the scaling functions are very similar; the main difference is that the scaling 
function decreases slightly more slowly for n = 5 than n = 4. This is reflected in the 
parameters of the fitting function for the small-x range: the amplitudes a and /3 tend to 
decrease as n increases (table 1). 

We turn to the discussion of the simulation results for one- dimensional systems. We shall 
describe the relevant behavior in Fourier space. Real-space data has been presented in ||. 
Our results show that in one-dimensional systems the asymptotic behavior of the structure 
factor also has a stretched exponential form, but the fitted exponents 5 are larger than those 
of the corresponding two-dimensional models, and close to unity for n = 4 and 5. 

In Fig. 5(a), we present the simulation results for the scaling function g(q) of the structure 
factor for the one-dimensional Heisenberg Model (n = 3). The continuous curve is the result 
of the analytical approach described in section IV. The analysis of the asymptotic behaviour 
gives an exponent 5 = 0.79 for the stretched exponential. Fig. 5(b) shows the curve of 
ln[g(g)] versus q 6 , where the linear behaviour is clearly observed. Similarly, we present the 
corresponding plots for the n = 4 model in Fig. 6, where the measured exponent is now 
5 = .98, while for n = 5 we obtain 5 = 1.02 as is shown in Fig. 7. The values of S for the 
last two models are so close that in practice it is difficult to distinguish between them. They 
are also close to the value unity obtained from the approximate large-n equation discussed 
in the following section. 

It is clear from the results for the one-dimensional models that the scaling function in 
real space is not a gaussian function, despite the good real-space fits to this form obtained 
in ||. Moreover, the (effective) exponents 5 for the n = 4, 5 models are bigger than for the 
corresponding models in two dimensions. Therefore, we have evidence that the exponent 5 
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increases with n, while it seems to decrease with d. Note that the gaussian result obtained 
for n = oo corresponds to 5 = 2, so the results presented here for the structure-factor tail are 
actually quite far from that limit. The analytical treatment presented in the next section 
gives some indication of why this might be expected. In particular, it suggests that the 
structure factor is dominated by a simple exponential for q — > oo at fixed large n, while the 
familiar gaussian form is recovered as n — > oo at fixed q. 

We conclude this section by discussing briefly some alternative fitting forms for the 
structure factor tail. First, however, we note that the stretched exponential form ( |3 . 1\) 
describes the tail well over at least 10 decades of S(k) in all cases. Of course, this represents 
a much smaller dynamic range (1 to l| decades) in the scaling variable q = kL(t). Motivated 
by the analytical result [equation ( }4.7|) below] of the approximate large-n theory, other fitting 
forms were tried for d = 2 (the agreement with the large-n theory already being good for 
n = 4 and 5 in d = 1). A direct fit of ( f4.7|) does not work well for d = 2. Allowing for a general 
power-law prefactor, g(q) = Aq~ x exp(-bq), gives a reasonable fit, but with very large values 
for x — 5.6 for n = 4, and 6.7 for n = 5. Fixing x = 1/2, but allowing for a general stretched 
exponent 5 again gives a reasonable fit (with 5 ~ 0.68 for n = 4 and 0.70 for n = 5), but 
over a significantly reduced range of q. For these reasons we prefer the unmodified stretched 
exponential ( [3.1|) as giving the simplest and most convincing description of the large-g data, 
at least for these small values of n in d — 2. Of course, it is quite possible that the form 
( |4.7|) will fit the data well at larger values of n. 

IV. ANALYTICAL TREATMENT 

In an attempt to gain some analytical insight into the structure factor asymptotics, we 
start from an approximate equation of motion for the pair correlation function derived using 
the gaussian auxiliary field approach pioneered by Mazenko ||20|| . We then make, for reasons 
that will become clear, the further simplification of retaining only the leading nonlinearity 
as n — > oo. The resulting equation is then finally used to extract the asymptotics of g(q). 



10 



The GAF method for vector fields has been discussed in some detail elsewhere. We refer 



the reader to the original papers and a recent review for a full exposition. The 
essence of the method is a mapping from the original field variable <ft to an 'auxiliary field' 
rh. The function 0(m) satisfies the equation (1/2) Y%=i d 2 (p/d^ rli = dV/dcj), where V{4>) is 
the potential in the Ginzburg-Landau function (4). With the boundary conditions 0(0) = 0, 
and 4>{m) — *> rh/\rh\ for \rh\ — > oo, this equation for 0(m) represents the equilibrium profile 
function for a spherically symmetric topological defect, with |m| representing distance from 
the defect. 

The (uncontrolled) approximation that m is a gaussian field, and the imposition of the 
scaling form (1), leads eventually to the self-consistent equation P1JT4 



r+ (^ + *) / + ^0 (4,) 

for the scaling function f(x), where primes indicate derivatives. In ( f4.1| ) 7 is the normalized 
correlator of the auxiliary field, 7 = (m(l) ■ m(2))/[(m 2 (l))(m 2 (2))] 1 / 2 , where '1' and '2' 
represent the space-time points x l5 t and x 2 ,t, and the function C{j) is given by 



CM = £ 



B 



n + 1 1 



2 



2 ' 2 

where B(x,y) is the beta function, and F(a,b; c; z) the hypergeometric function. The con 



stant A in ( |4.1|) has to be adjusted so that f(x) vanishes sufficiently fast at infinity ||14|| . 

As should be clear from the above discussion, (|4.1|) only really makes sense for n < d, 
based as it is on the presence of singular topological objects whose positions are defined by 
the zeros of the field or, equivalently, by the zeros of rh. Indeed, the function C(j) has 
inbuilt structure that generates the Porod tail associated with such defects. Specifically, in 
the short-distance limit, where 7 — > 1, the hypergeometric function in ( |4.2| ) has a singular 
contribution of order (1 — 7 2 ) n / 2 (with a logarithmic correction for even n). Since 1 — 7 2 ~ x 2 
for small scaling variable x, this singular term is of order x n (again, with a logarithm 
for even n), leading to the power-law tail g(q) ~ q-( d + n ) i n Fourier space. Within the 
GAF approach, this tail is obtained for all n and d. For n > d + 1, however, neither 
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singular topological objects nor nonsingular topological textures exist, so the GAF result is 
qualitatively incorrect. Indeed, this is to be expected since the GAF approach is specifically 
designed to build in the defect structure. 

So what should one do when there are no defects? We have seen that the usual GAF 
approach always give a Porod tail, for any n and d: this is unphysical for n > d + 1, since 
the tail is a consequence of the presence of topological defects. One way around this impasse 
is to artificially approximate the full GAF equation Q4.1]) by the form valid for n — > oo. In 
this limit jdC '/dj = f + f 3 /n + 0(l/n 2 ), and Q4.1|) becomes, correct to 0(l/n), 



r+ (^ + £),^( /+ I /3 ).0. (4.3) 

This step, admittedly ad-hoc, has the desired effect of eliminating the unwanted (for 
n > d + 1) short- distance singularity in f(x). Eq. ( f4.3|) is the nonconserved version of 
the equation introduced by BH to study the crossover from multiscaling to simple scaling in 



the asymptotic dynamics of a conserved vector field at large but fixed n ||12|| . Both conserved 
and nonconserved versions have recently been studied numerically [ PH . 

To extract analytically the large-g behavior, we perform a (ci-dimensional) Fourier trans- 
form of (|4.3|). The resulting equation for g(q) = J d d x f(x) exp(iq • x) is 



'\ + q^jg{q) + \g'{q) = \ <J}(q) + B(q)) 



(4.4) 



where 



B(q) = -( d d x f{x) exp(*q • x) . (4.5) 
n J 



If we assume an asymptotic form g(q) ~ q v exp(—bq ' ), with S < 2, then ( f4.4|) gives 

B(q) ^ ^-g{q) , q^oo. (4.6) 



In the Appendix, we show that consistency with (4.6) requires 5 = 1 and v = (1 — d)/2, 



i.e. 



g(q) -> Aq (1 - d)/2 exp(-bq) , q -> oo . (4.7) 
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In real-space this implies that the function f(z) has simple poles in the complex z plane at 
z = ±ib. The value of b is not determined by this argument; instead one can derive (see 
Appendix) the relationship 



A = (lQn z n/X) (27tb) a - 1 (4.8) 

between b and the prefactor A in the asymptotic form ( |4.7| ). The existence of these simple 
poles in real space also follows directly from the real-space equation Q4.3|) . If one assumes a 
singularity of the form (z — Zq)' 1 , with 7 > 0, then balancing the dominant terms f" and 
Xf 3 /2n gives immediately 7 = 1, i.e. a simple pole. The position zq is not determined, but 
the residue C of the pole is given by C = =p(4n/ A) 1 / 2 , where the two values correspond 
to the poles z = ±ib. Using this result for C, one can readily recover ( f4.8| ) by contour 
methods, e.g. for d — 1 one has 



g(q) = / dx f (x) exp(iqx) 

J -co 

-> 27r(4n/A) 1/2 exp(-6g) , q -> 00 , (4.9) 

where the second line, equivalent to ( f4.8| ) for d = 1, was obtained by closing the contour in 
the upper half-plane. 

The approach outlined above gives the relation fl4.8|) between A and b, but does not 
determine b explicitly. We now give a heuristic argument that b ~ (Inn) 1 / 2 for large n. First 
we make an observation concerning the value of A. Equation ( f4.4|) with q = gives 

a = i m 



2^(0) +5(0) 



d 



1 



1 + - / d d xf-\x)/ / d d xf(x) 



-1 



(4.10) 



2 

In particular, A = d/2 for n = 00. For n = 00, therefore, ( |4.4| ) becomes q 2 g(q) + (q/4=)g'(q) = 
0, with solution (7(g) = (S^)^ 2 exp(— 2g 2 ) [the prefactor being fixed by the condition /(0) = 
1]. For n large but finite, on the other hand, we have seen that the asymptotic form 
is g(q) ~ exp(— bq). The crossover between these two forms presumably occurs at some 
q = q*(n), with g(q) ~ exp(— 2q 2 ) for 1 q q*, and g(q) ~ exp(— bq) for g ^> q*. 

13 



Matching these two forms at q = q* gives q* ~ b. Next we evaluate B(q) in the region q <C q*. 
Here q(q) ~ (87r) d / 2 exp(— 2g 2 ), so /(s) ~ exp(— x 2 /8), giving 5(g) ~ (1/n) exp(— 2g 2 /3). 
However, this decays more slowly with q than the other terms in ( |4.4j ), which fall off as 
exp(— 2q 2 ). So the term involving B(q) (evaluated for q <C q*) becomes comparable with the 
other terms when (1/n) exp(— 2g 2 /3) ~ exp(— 2q 2 ), i.e. when q ~ (Inn) 1 / 2 . This suggests 
q* ~ (lnn) 1//2 , and therefore b ~ (Inn) 1 / 2 . The numerical data of Castellano and Zannetti 
|l~3| certainly show that b increases extremely slowly with n. 

To compare this approximate theory with our simulation data, we have solved Ql.l|) 



numerically for d— 1, n— 3, 4, 5, using the procedure described in JTj]. The Fourier trans- 
form was then taken numerically, and the 'best fit by eye' to the structure-factor data was 
obtained by adjusting the timescale in the theoretical curves, giving the results shown by 
the continuous curves in Figures 5(a), 6(a) and 7(a). The corresponding log-linear plots, 
which reveal the large-g behaviour more clearly, are shown in Figures 8(a)-(c). As might be 
expected, equation ( |4.1| ) [or its Fourier transform (|4.4|)1 does not describe the data quantita- 
tively over the whole range of q = kL(t), but it does give a qualitatively correct description. 
There is an early parabolic region, corresponding to a gaussian form for g(q), which then 
gives way to a slower decay that, at least for n = 4 and 5, is consistent with the simple 
exponential form predicted by ( [4.4|) but with a different coefficient b in the exponent. Given 
that the theory is, at best, a large-n theory we regard these results as encouraging. The 
n = 3 data, however, and the d = 2 data, do not seem to fit a simple exponential, at least 
for the range of q that we have been able to explore. (This is of course implicit in the values 
5 < 1 obtained for these systems from Figures 2,4 and 5.) It may well be that considerably 
larger values of n are needed in d = 2 than in d = 1 for the large-n asymptotics to become 
apparent. 

The above derivation of an exponential tail was specific to nonconserved fields. What can 
we say for conserved fields? The fundamental equation of motion for this case is obtained 
from the TDGL equation ( |2.1| ) by the replacement T — * —TV 2 . Applying the GAF method 
to this equation, imposing the scaling form (|1.1|) [but with L{t) = t 1//4 for conserved fields], 

14 



and taking the Fourier transform, leads to [12 



| + /) 9(q) + ^'(q) = ^ (*(?) + B(q)) , (4.11) 



instead of ( f4.4j ). [The definition ( |4.5| ) of B(q) differs by a constant from that used in [|T2 
where A was written as 2g^, g m being the position of the maximum of g(q) for large n.] 
Assuming the asymptotic form q(q) ~ q v exp(— bq s ) for q — > oo, ( |4.11| ) gives ( |4.6|) once more, 
provided 5 < 4. Then our previous arguments apply, and the asymptotic form (|4.7| ), with 
A and related by (|4.8|) , are recovered. This approach therefore predicts that the structure 
factors for conserved and nonconserved systems will have the same asymptotic forms, at 
least within the context of the BH truncation. The same conclusion was drawn from recent 



numerical solutions of the BH equation [13 . 



V. CONCLUSION 

In summary, we have studied the dynamics of phase ordering for models without topo- 
logical defects in one and two dimensions. We find that scaling is achieved with the growth 
law Lit) = t x l 2 . The tail in the structure factor is well fitted by a stretched exponential form. 
For the two-dimensional systems, table 1 summarizes the relevant parameters describing the 
fits in real and Fourier space. In contrast to systems with singular defects (n < d), where 
the generalized Porod form g(q) ~ q-( d + n ) for the structure factor tail is a consequence of 
the defect structure, and is independent of the presence or absence of conservation laws, in 
systems without defects the functional form does, apparently, differ for conserved and non- 
conserved systems. We have shown, for example, that for the particular case of the n = 4 
model in two dimensions the tail is well described, over the range of q accessible to us, by a 
stretched exponential with exponent 5 = .435, differing from the result for the same model 
with conservation studied by RC ||10|| , who found 5 ~ 1.7. 

Within the 'toy' equation of Bray and Humayun fll^ |, however, we have shown that 



the true asymptotics are the same for conserved and nonconserved dynamics. Of course, 
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the BH equation is at best a large-n theory, and the numerical results for nonconserved 
and conserved dynamics may converge as n is increased. A related question is whether 



the exponents 5 measured here and in jTD| are genuine asymptotic exponents, or effective 
exponents whose values will change as the range of q over which the fit is made is moved 
to larger q. More extensive simulations may cast some light on this issue. The 'universal' 
(independent of conservation laws) Porod tail behavior obtained for n < d is geometrical in 
origin, being a consequence of the field structure induced by singular topological defects [0 . 
As yet, however, we have no corresponding physical picture in the absence of topological 
defects. 

It is interesting that, within the simple model of equation ( |4.1|) , the exponent 5 jumps 
discontinuously from 5 = 2atn = ooto5 = l for n large but finite. More precisely, one can 
say that 5 = 2 corresponds to the limit n — > oo at fixed, large q, while 5 = 1 corresponds to 
q — > oo at fixed, large n. We have argued that the crossover between these limiting forms for 
fixed, large n occurs at q ~ (lnn) 1//2 . This change of behavior depending on the order of the 
limits is reminiscent of the result obtained from the conserved version of ( [lip , where a novel 



'multiscaling' behavior is obtained for n — > oo at fixed, large t [IT], while simple scaling is 



recovered for t — > oo at fixed, large n WM. For the nonconserved case, one always has simple 



scaling. For both conserved and nonconserved fields, however, the asymptotics of g(q) are 
sensitive to whether n is large or truly infinite. This rules out, for example, exploring the 
asymptotics by expanding around the large-n solution in powers of 1/n. 
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VI. APPENDIX 



In this Appendix we use ( |4.6|) to derive the asymptotic form (|4.7|) for g(q). From the 
definition ( |4.5|) we have 

B(q) = - / / g(p)gQt)g{q - p - k) , (6.1) 
where J p = J d d p/(27c) d . Inserting the asymptotic form 

g(q)^Aq l/ exp(-bq s ) , (6.2) 

gives 

B(q) - — / / F(p, k, q) exp[-&£(p, k, q)] , (6.3) 
n Jp Jk 

where 

F( P ,k,q) = | P nkriq-p-kr 

E(p, k, q) = |p| 5 + |k| 5 + |q - p - k| 5 . (6.4) 

We now scale out the g-dependence through the changes of variable p = qu, k = qv, q = qe, 
where e is a unit vector. Then 

B(q) = — q 2d + 3 » f f F(u,v,e) exp[-bqE(u, v, e)] . (6.5) 

For q — > oo, we can attempt to evaluate the u and v integrals using the method of 
steepest descents. This requires minimizing the function E(u, v, e). The points requiring 
consideration are the symmetry point, u = v = e/3, and the points u = = v and two 
similar points obtained by permuting u, v and e — u — v. The corresponding values of E 
are E(e/3, e/3, e/3) = 3 1 -" 5 , and £(0,0, e) = E(0,e,0) = E{e,0,0) = 1. Thus for 5 > 1, 
the symmetry point minimizes E, giving B(q) ~ exp(— 3 1 ~ s bq). But this form violates 
the asymptotic relation (|4.6|) , according to which B(q) and g(q) must decay with the same 
exponential factor, so 5 > 1 is ruled out. 
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For 5 < 1, the smallest E is unity, obtained when two of u, v, and e — u — v vanish. 
So this case is apparently consistent with ( f4.6| ). However, the integral is now dominated 
by points where two of the momenta p, k, and q — p — k vanish. This invalidates the use 
of the asymptotic form for g(q) in the evaluation of B(q), so the derivation of a stretched 
exponential form is not internally consistent for 5 < 1. 

This leaves 5 = 1. For this case all points of the form u = ae, v = fie, with < a < 1 
and < fi < 1 — a, give E = 1, so one has to integrate over all such points. Writing 
u = ae + u_i_, v = fie + vj_, expanding E to quadratic order in u^, vj_, and carrying out 
the integrals over u^, vj_, gives after some algebra 



B{q) = (A 3 /4ir 2 n) q d+1+3u exp(-bq) I(d, v)/(2nb) 
I(d, v) = f 1 da [~ a dp [ap{l -a- fi)Y+(d-i)/2 _ 



jo Jo 



d-l 



(6.6) 



But ( f4.6|) implies, asymptotically, 



B(q) = (2A/X) q 2+u exp(-bq) . 



(6.7) 



Comparing ( p.dp and ( |6.7| ) gives v = (1 — d)/2 and Eq. ( ^4.8| ) for the amplitude A. 
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FIGURES 



FIG. 1. (a) Real-space correlation function C(r) for a 2D system with n = 4 components 
as a function of distance r for several times t. The data were obtained from lattices of 
size 256 x 256, averaged over 20 different initial conditions, (b) Demonstration of dynamic 
scaling, C(r,t) = f(r/L(t)), with L(t) = t 1 ' 2 . 

FIG. 2. (a) Scaling structure factor g(q) = [L(t)]~ 2 S(k, t) as a function of scaled mo- 
mentum q = kL(t) for a 2D system with n = 4 components for lattices of size 256 x 256 
(averaged over 20 different initial conditions), (b) Demonstration of the stretched exponen- 
tial behaviour, plotting In g(q) against q 5 with 5 = 0.435. The line is included as a guide to 
the eye. 

FIG. 3. Same as Fig. 1 but for the n = 5 model. 

FIG. 4. (a) Same as Fig. 2 but for the n = 5 model, (b) The power 5 = 0.613 was found 
to give the best linear relation between In g(q) and q s . 

FIG. 5. (a) Scaling structure factor g(q) = [L(£)] -1 S(k, t) as a function of scaled momen- 
tum q = kL(t) for a ID Heisenberg model [0(3) model] for lattices of size 16384 (averaged 
over 100 different initial conditions). Continuous curve: result of the approximate analytical 
treatment described in the text, (b) Demonstration of the stretched exponential behaviour, 
plotting \ng(q) against q 5 with 5 = 0.79. The line is a guide to the eye. 

FIG. 6. (a) Same as Fig. 5, but for the n = 4 model, (b) The power 5 = 0.98 is found 
for this model to give the best linear relation between lng(g) and q 5 . 

FIG. 7. (a) Same as Fig. 5 but for the n = 5 model, (b) Same as Fig. 5(b), but with 
5 = 1.02. 

FIG. 8. Log-linear plot of the scaled structure factor against scaled momentum for the 
ID systems: (a) n = 3 (b) n = 4 (c) n — 5. In each case the continuous curve is the result 
of the approximate theory of section IV. 
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TABLES 

TABLE I. Parameters determined from fitting of the simulation results for the two-dimensional 
systems : the small-x solution (x < 0.5) in the scaling function with the form f(x) = 1 — ax 2 + /3x 4 
and the asymptotic analysis of scaling structure factor in term of the stretched exponential 
g(q) ~ exp(-bq s ) 
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Tail 
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5 
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0.435 




8.19 
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0.8417 


0.613 




4.39 
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